NASA Contractor Report 181913 
ICASE Report No. 89-65 


ICASE 


ON THE COEFFICIENTS OF DIFFERENTIATED 
EXPANSIONS OF ULTRASPHERICAL POLYNOMIALS 


- Andreas Karageorghis 
Timothy N. Phillips 


Contract No. NAS 1-1 8605 
September 1989 


Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 

Operated by the Universities Space Research Association 

(NA^A-CR-1 81913) ON THE COEFFICIENTS OE N<?0~i0o37 

D I FF E R t NT I AT t D FXPANSTONS OF ULTRASPHFRICAL 
POLYNOMIAL S Final Report (TCASF) Up 

CSCL 12 A Unci as 

03/64 0234436 

l\!/\S/V 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665-5225 



ON THE COEFFICIENTS OF DIFFERENTIATED 
EXPANSIONS OF ULTRASPHERICAL POLYNOMIALS 


Andreas Karageorghis 1 
Mathematics Department 
Southern Methodist University 
Dallas, TX 75275-0156 

and 

Timothy N. Phillips 2 
Department of Mathematics 
University College of Wales 
Aberystwyth, SY23 3BZ, U.K. 


ABSTRACT 

A formula expressing the coefficients of an expansion of ultraspherical polynomials which 
has been differentiated an arbitrary number of times in terms of the coefficients of the original 
expansion is proved. The particular examples of Chebyshev and Legendre polynomials are 
considered. 
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1. Introduction 


The current popularity and interest in spectral methods is mainly due to their superior 
approximation properties compared with other methods of discretization. One can show that 
if the eigenfunctions of a singular Sturm- Liouville problem are used as basis functions then 
the rate of decay of the expansion coefficients is determined by the smoothness of the function 
being expanded and not by any special boundary conditions satisfied by the function. If the 
functions of interest are infinitely differentiable then the nth coefficient a n decreases faster 
than any finite power of £ as n — » oo [ 2 ]. 

Spectral methods have been extensively used in the solution of boundary value problems 
and in computational fluid dynamics [1, 2]. In most of these applications use is made of for- 
mulae relating the expansion coefficients of derivatives appearing in the differential equation 
to those of the function itself. In fact, the coefficients of successive derivatives of a function 
are related by a recurrence relation which can be solved numerically in an efficient and sta- 
ble way. This is a distinguishing feature of the pseudospectral method [ 1 ]. For the Galerkin 
and tau variants of the spectral method explicit expressions for the expansion coefficients of 
derivatives in terms of the expansion coefficients of the solution are required. In this paper, 
we obtain such expressions when the basis functions are ultraspherical polynomials. This 
work generalizes earlier results for the Chebyshev and Legendre polynomials [3, 4] and is the 
converse of a result on the coefficients of integrated expansions of ultraspherical polynomials 


[5]- 

In Section 2 we give relevant properties of ultraspherical polynomials and in Section 3 
we describe how they are used to solve boundary value problems with the tau method. In 
Section 4 we prove the main result of the paper which is an expression for the coefficients 
of a general order derivative of an expansion in ultraspherical polynomials in terms of the 
coefficients of the original expansion. Concluding remarks are given in Section 5. 


2. Properties of Ultraspherical Polynomials 

The ultraspherical or Gegenbauer (see [ 6 ]) polynomials {Pj[ 7 )(x)}^l 0 are orthogonal on 
the interval [- 1 , 1 ] with respect to the weight function Wy(x) = (1 — x 2 ) 7 , i.e. , 

j\ 1 - x 2 yp£\x)pW{x)dx = 0, m^n. (1) 

In order that the weight function is non-negative and integrable we assume 7 > — 1. The 
ultraspherical polynomials are eigenfunctions of the following singular Sturm- Liouville prob- 
lem 

(1 — x 2 )<f>"(x) — 2(7 + 1 )x<f>'(x) + n(n + 27 + l)<^>(x) = 0 . ( 2 ) 

A consequence of this is that spectral accuracy can be achieved for expansions in ultraspheri- 
cal polynomials. This class of orthogonal polynomials is a subclass of the Jacobi polynomials. 

Let /(x) be an infinitely differentiable function defined on [-1,1]. We may represent /(x) 
in the form 

OO 

/(*) = £ <■»#'(*). ( 3 ) 

71 — 0 
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where Pj?\ x) is the ultraspherical polynomial of degree n. Further, let a£0 denote the 
ultraspherical coefficients of the qth derivative of /(x) then it can be shown [5] that the 
coefficients of successive derivatives of /(x) satisfy the recurrence relation 

0 («-i) = (^ + 27) (,) (n + 7 + 1) (,) 

n (n + 7 )(2n + 2 7 - 1) n " 1 (» + 2 7 + l)(2n + 2 7 + 3) n+1 ‘ 

This recurrence relation may be solved to yield 


= (2n + 2 7 + 1) 


r(n + 2 7 + 1) 

T(n + 7 + l) 


E 


k=l 


r(n + 7 + 2 k) (,_ 1 ) 

T(n + 2 7 + 2fc) 0n+2fe-1 ’ 


( 5 ) 


for q > 1 where aj, 0 ' = a n and T(-) is the Gamma function. 


3. The Tau Method for a Model Problem 


Consider the boundary value problem 

u"(x) + Au(x) = g(x), xe[— 1, 1], (6) 

u(±l) = 0, (7) 

where A is a real scalar. Suppose that we approximate u(x) by a truncated expansion of 
ultraspherical polynomials 

Un(x) = X) “nPnK*)- ( 8 ) 

n= 0 

We seek to determine a n (0 < n < N) using the tau method. 

The coefficients a n are chosen so that ujv(x) satisfies 

<(x) + \u N (x) = g N (x), (9) 

where 

9n(x) = ( 10 ) 

n=0 

Substituting (8) and (10) into (9), multiplying by w-,(x)Pffl(x), 0 < m < N — 2, and 
integrating over [-1,1] yields 

a( m + Aa m = d m) 0 < m < N — 2. (11) 

Two further equations are obtained from the boundary conditions by inserting (8) into (7): 


N 


E*n(±l)’ 


r(n + 7 + 1 ) 
r(7 + 1) 


= 0 . 


( 12 ) 


Thus (11) and (12) comprise N + 1 equations for the N + 1 unknown coefficients a n (0 < n < 

N). 
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The coefficients, of the second derivative of the approximation Ufj(x) are related to 
the coefficients, a^, of ujv(x) by invoking (4) with q = 1 and q = 2. In the next section 
we show how the coefficients of any derivative may be expressed in terms of the original 
expansion coefficients. This allows us, for example, to replace in (11) by an explicit 
expression in terms of the a n . In this way we can set up a linear system for a n (0 < n < N ) 
which may be solved using standard direct solvers. 

4. The Theorem and its Proof 
Theorem 


Let f(x ) be an infinitely differentiable function defined on the interval [-1,1]. The coef- 
ficients oj*) of an expansion of ultraspherical polynomials of the qth derivative of /(x) are 
related to the coefficients o n of an expansion of ultraspherical polynomials of f(x) by 



(2n + 2 7 + l)r(n + 2 7 + l)y, T(n + 7 + 2k + q - 1) (k + q - 2)! 
2 q ~ 2 (q — l)!r(n + 7 + 1) ~r(n + 2 7 + 2* + g-l) (*-1)! 


r(2n -j- 2 7 -|- 2k -|- 2 q — 2) 
r(2n -f- 2 7 -f- Ik -(- 1) 

for all n > 0, q > 1. 


T(n + 7 + k + 1) 
r(n + 7 + k + q - 1) n+2k+q 2 


(13) 


Proof 


For q = 1 (13) reduces to (5) when q = 1. We proceed by induction and assume that the 
theorem is valid for q. We need to show that it is valid for q + 1, i.e., 

f0+n (2n + 2 7 + l)T(n + 2 7 + 1) “ T(n + 7 + 2k + q)(k + q - 1)T 

“ 2«- 1 g!r(n + 7 +l) ^ T(n + 2 7 + 2fc + g) (u) 

(2 n + 2 7 + 2k + 2g)r(n + 7 + k + 1) 

(J k - l)!r(2n + 2 7 + 2ifc + l)r(n + 7 + k + “) a "+ 2fc +’- 1 ’ 


We replace q by q + 1 in (5) and then substitute for using (13), which is the inductive 
hypothesis, to give 


a£ +1 > = (2n + 2 7 + l) 


r(n + 2 7 + 1) 
T(n + 7 + 1) 


E 


fc=i 


r(n + 7 + 2fb) [(2n + 2 7 + 4A:-l)r(n + 2 7 + 2^) 
r(n + 2 7 + 2Jfc) [ 2 q ~ 2 {q - 1)! T(n + 7 + k ) 


” r(n + 2 ± 2fc + 2j + q - 2 )(j +q- 2)! 

X ,tl r(n + 2 7 + 2A: + 2i + 9-2)(j-l)! 

r(2n + 2 7 + 4 k + 2 j + 2 q- 4)f(n + 7 + 2k + j) 
r(2n + 2 7 + 4Jfc + 2j - l)T(n + 7 + 2 k + j + q- 2)‘ ln+2k+2i+ q ~* 

(2n + 2 7 + l)r(n + 2 7 + l) ” f* Tjn + 7 + 2k + 2; + q - 2) 

2 q ~ 2 (q — l)!T(n + 7 + 1) ^ + *7 + ™ J k F ( n + 2 7 + 2k + 2j + q - 2) 
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Q + <Z — 2 )!r( 2 n -j- 27 + 4A: + 2 j + 2 ^ — 4) T(n + 7 + 2 k + j) 

(j - l)!r(2n + 7 + 4 k~+ 2 j - 1) X r(n + 7 + 2fc + j + q - 2 ) 0n + 1 2 *+ 2 >+*- 3 ' 

If we write A + j = l + 1 then we may write in the form 




( 2 n + 27 + l)r(n + 27-4-1) ^ T(n + 7 + 21 + q) 
2 q ~ 2 (q — l)!r(n + 7 + I) ^ T(n + 27 + 21 + q) 


r 1 


XX 2 n + 27 + 4A - 1) 

Jb=l 


x 


(£-k + q-l)\ 


r( 2 n + 27 + 2k + 21 + 2q - 2 )r(n + 7 + & + I + I) 
r(2n -(- 27 + 2k + 21 4 - l)r(n + 7 4” ^ ~ h 9 — I) 


°n+2/+g-l- 


In order to complete the proof of the theorem we need to prove the following lemma: 


(15) 


Lemma 


T(2n + 2 , + 4 * _ n (*-fc+t-l) !r ( 2 ™ + 27 + 2 < ; + 2 l + 2 <|- 2 )r(n+ T + *+<+l) 
k = i {£ — &)!r( 2 n + 27 + 2k + 2t + l)r(n ■j-'y-hk-i-£-i-q — 1 ) 

= 1 (^ + g-l)!r(2n + 27 + 2/ + 2g)r(n + 7 + ^+l) 

2 q (l - l)!r( 2 n + 2 7 + 21 + l)T(n + 7 + t + q) ' 

for all values ofn>0,£>l,g>l. 


Proof 

When 1 = 1 both sides of (16) are equal, taking the value 

{q — l)!r(2n -+■ 27 + 2 q + 2)T(n + 7 + 2) 

2 r(n + 7 + q + l)r( 2 n + 2 7 + 3 ) ’ 

and in this case (16) holds for all values of n > 0,g > 1. The proof of the lemma proceeds 
by induction on l. We need to prove that (16) is true when l is replaced by l + 1, i.e., 

vV n 1 | . ■ ,' (*-* + 9)!r(2n + 2 7 + 2 k + U + 2g)r(n + j + k + l + 2) 

S * + l)!r(2n + 2 7 + 2it + 2< + 3)r(n + 7 + *: + < + ?) 

( 17 ) 

_ 1 (^ + 9)! r(2n + 27 + 2 q + 2)r(n + 7 + ^ + 2) 

2 q l\ r(2n -(- 27 ~h 2£ -f- 3 )r(n + 'y + £ + q 1) 

If we replace n by n + 2 in (17) we obtain 

y'/ 2 n 1 2-7 4- + 3*\ ~ ^ 4 £ ~ r(2n + 27 + 2& + 21 + 2 q + 2)r(n + 7 -I-A-M + 3) 

££ 7 ' (/-A)! r( 2 n + 2 7 + 2 A + 2 ^ + 5 )r(n + 7 + k + i + q + 2) 

1 {l + q - 1)! r(2n + 2 7 + 21 + 2 q + 4 )r(n + 7 + ^ + 3 ) 

2 q (£ — 1)! r(2n + 27 + 2/ + 5 jr(n + 7 + / + g + 2)* 

(18) 
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If we set p — k + 1 then (18) becomes 


(l-p + <z)!r( 2 n + 2 7 + 2 p + 2l + 2 g )r(n + 7 +p + l + 2) 
2^{2n + 2 7 + p > (i_ p + l)!r(2n + 2 7 + 2p + 21 + 3)T(n + 'y+p + l + q ) 

_ 1 (£+g_- l)!r(2n + 2 7 + 21 + 2 q + 4)f(w + 7 + £ + 3) 

~2g (/- l)!r(2n + 2 7 + 2^ + 5)r(n + 7 + ^ + 9 + 2) 

Now the left-hand side of (17) may be written as 

^ + g -l)!r( 2 n + 2 7 + 2^ + 2 g + 2)r(n + 7 + ^ + 3) 

[ + 7 + j £!T(2n + 27 + 21 + 5)r(n + 7 + £ + g + l) 


i , (/-fe + g)!r(2n + 27 + 2fc + 2/ + 2 g )r(n + 7 + < ! + < + 2) 

+ hi + 27 + ’ {l - k + l)!r(2n + 27 + 2k + 21 + 3)r(n + 7 + fc + * + g)’ 

which, in view of (19), is equal to the following 

, „, (£ + g-l)!r(2n + 2 7 + 2^ + 2 g + 2)r(n + 7 + ^ + 3) 

[ 2 n+ 7+ ) ^!r(2n + 27 + 2^ + 5)r(n + 7 -M + 9 + 1) 

1 (l + q- l)!r(2w + 27 + 21 + 2q + 4)r(n + 7 + 1 + 3) 

+ 2q (l-l)\T(2n + 2'y + 2l + S)r(n + 'y + l + q + 2) 

1 (l + q - l)!r(2w + 27 + 21 + 2q + 2)T(n + 7 + 1 + 3) 

“ 2q {l- l)!r( 2 n + 27 + 2 £ + 5)r(n + 7 -M + g+ 1 ) 

(2 n + 2 7 + 3)2g (2 n + 2j + 2i + 2g + 3)(2n + 2 7 + 21 + 2q + 2) ' 

X t + (n + 7 + ^ + 9 + 1 ) 

1 {l + q- l)!r(2n + 27 + 21 + 2q + 2 )r(n + 7 + i + 3) \ 2 ( 2 n + 2 7 + 21 + 3)(g + O' 

~Tq [t - l)!r(2n + 2 7 + 2i + 5)r(n + 7 + t + q + 1) [ l 

1 (l + g)! r( 2 n + 27 + 21 + 2q + 2 )r(n + 7 + l + 2 ) 

" 2q l\ r( 2 n + 27 + 2 ^ + 3)r(n + 7 + ^ + 9 + l) ’ 

which is the right-hand side of (17). This completes the inductive step and proves the lemma. 

Application of the lemma to the expression in square brackets on the right-hand side of 
(15) gives (14) which completes the proof of the theorem. O 

The particular expressions for the Chebyshev and Legendre polynomials may be derived 
as special cases of the theorem. We give these as corollaries to the main theorem. 


Corollary 1 ([4]) 

If /(x) is an infinitely differentiate function on [-1,1] then the Legendre coefficients, 
of the qth derivative of f(x ) are related to the Legendre coefficients, a n , of f(x) by 



(2n + 1) A (k + q - 2)! (2n + 2k + 2q- 3)!(n + *)! 

29 - 2 ( g _ 1)!^ (jfc-1)! (2n + 2k)\{n + k + q-2)\ an+2k+q - 2 ’ 


( 20 ) 
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for all n > 0, q > 1. 

Proof 

Simply set 7 = 0 in (13). 

Corollary 2 ([3]) 

The corresponding Chebyshev derivative coefficients are related by 

,(o) 2 g ^2, [k -\- q — 2)!(n + k + q — 2)!. 

^ n (q — 1)! (Jfc- l)!(n + Jfe- 1)! ^ + 2k + q ~ 2 ) b n+2k +q -2, (21) 

where Co = 2, c* = l(n > 0). 

Proof 

We consider the case n > 0. The ultraspherical polynomials with 7 = — | are related to 
the Chebyshev polynomials by 


Pi-' m (x) = 2~ ln ( 2 ") r„(x), 


n > 1, 


The ultraspherical coefficients, of the qth derivative of f(x) when 7 = — ^ are therefore 
related to the b$ by 


4 ,) = for «u ? > 0. 

(2n)! 

Setting 7 = — | and substituting from (13) into (22) we obtain 

2nr(n) 


( 22 ) 


2 = 


1 “ (k + q-2)\T{n + 2k + q- f) 

i\ 2s 


(2 n)\ n 2 g ~ 2 (q - 1)! r(n + I) £ (* - l)!r(n + 2* + q - 2) 


r( 2 n + 2k + 2q- 3 )r(n + * + 5) n2n+4k+2q _ 4 [(n + 2k + q- 2 )!] 2 t 

^ T \ -r-, / , Qv * & !” I H TTT 0 , 


T(2n + 2&)r(n + A; + <7 — «) 


(2n + 4A: + 2q — 4)! 


^n+2fc+g— 2* 


After a little manipulation we may write this in the form 

bl q) = 


2«" 1 

(2 n)! 1 


r(n + 2 k + q — |)(n + 2 k + q — 2 )!' 

(9-1)! 

, n!r ( n + f). 


(2n + 4* + 2q - 4)! 


T(n + k + |) 

(2 n + 2k + 2q — 4)! 

r(2n + 2k) 

r(n + k + q - f ) 


( 23 ) 


>4fci 


Now application of the formula 


(2m)! 


2 2m rn! 


r(m + i) r(i) 
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to each of the terms in square brackets on the right-hand side of (23) yields (21) for n > 0. 
The special case n — 0 may be derived in a similar way. 

5. Conclusions 

This paper proves formulae associated with the differentiation of expansions of ultras- 
pherical polynomials and describes how they can be used to solve two-point boundary value 
problems. Formulae relating the coefficients in expansions which have been differentiated an 
arbitrary number of times to those in the original expansion are proved. 
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